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ABSTRACT 

We present new phase- resolved H and if-band spectroscopy of the ultra-short 
period magnetic cataclysmic variable EF Eri in its current, prolonged "low" state 
obtained using NIRI on Gemini-North, and NIRSPEC on Keck II. These new 
data show that the if-band spectrum of EF Eri appears to be dominated by 
cyclotron emission during the entire orbital cycle. The if-band spectrum of EF 
Eri is likewise dominated by cyclotron emission during most of an orbital period, 
but near binary phase 0.0, the secondary star spectrum may be visible. The lack 
of strong CO or CH4 absoprtion features, and the weakness of the water vapor 
features in this spectrum, however, suggests the possibility of peculiar abundances 
for Carbon and/or Oxygen. We have used the PHOENIX stellar atmosphere code 
to produce model brown dwarf spectra with non-solar abundances of Carbon, 
Nitrogen and Oxygen and achieved limited success in fitting the observed spectra. 
We conclude that strong, and highly variable cyclotron emission is responsible 
for the photometric variation previously reported for EF Eri. The nature of 
this cyclotron emission is complex: the if-band spectra show that the dominant 
cyclotron harmonic at phase 0.5 peaks at 1.65 //m, but at phase 0.0, the harmonic 
peaks near 1.72 /im. At phase 0.5, there is another cyclotron feature present that 
peaks in between the if and if bands (near 1.93 /^m), but at phase 0.0, no such 
feature is present. These data suggest that cyclotron emission from both poles is 
occurring. In the high state, the cyclotron emission has been modeled as coming 
from the pole that is oriented towards the secondary star. One interpretation for 
the phase 0.5 cyclotron emission is that it originates from the opposite pole. In 
its current ultra-low state, EF Eri reveals no outward signs of accretion (such as 
H I emission), but continues to have a few, strong cyclotron features. Thus, EF 
Eri joins the small group of magnetic cataclysmic variables whose accretion rate 



-4- 



is so low that they are in the "bombardment scenario" regime. 

Subject headings: binaries: close — stars: individual (EF Eridani) - stars: low-mass, 
brown dwarfs — stars: magnetic fields — stars: variables: general 
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1. Introduction 

Polars are magnetic cataclysmic variables that have white dwarf primaries with strong 
magnetic fields that dramatically effect the accretion fiow from their cool Roche lobe 
filling companions. Instead of forming a disk, the accreted matter fiows in a stream from 
near the Li point to one, or both, of the magnetic poles of the white dwarf primary (see 
Wickramasinghe & Ferrario 2000 for a complete review). For reasons not yet identified, 
the normal "high state" mass transfer rate (f» lO"-*^-^ Mq yr~^) ceases, and the system 
transitions to a "low state" . For the most extreme examples, accretion rates as low as 10~^^ 
Mq yr~^, or lower, have been inferred (Reimers et al. 1999, Schwope et al. 1999, Szkody 
et al. 2003). When in these low states, polars should present an excellent opportunity to 
examine the very low mass secondary stars of the shortest period cataclysmic variables 
(CVs). 

EF Eri is an ultra-short period (Porb = 81 min) polar that entered a low state in 
1997. Harrison et al. (2003) presented simultaneous multi-wavelength optical and infrared 
photometry for EF Eri in its low state in which they beheved they had detected the 
irradiated brown dwarf-like secondary star of this system. Their data showed that in the 
BVRIJ bandpasses, the light curves were dominated by emission from the white dwarf and 
some residual cyclotron emission, while in the H- and if-bands, the light curves exhibited 
large amphtude (AK = 0.8 mag) variations that repeated once per orbital cycle. Harrison 
et al. were able to model these multi- wavelength data with an extremely cool (Teff < 1,000 

^Partially based on observations obtained at La Silla (ESO proposal 69.D-0142(A)) 
^Visiting Astronomer, Kitt Peak National Observatory, National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foundation. 
^Visiting Astronomer, NTT Telescope at La Silla Operated by ESO 
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K) secondary star that is being irradiated by the primary. In this model, the irradiated 
hemisphere of the secondary star is heated to Teg 1,500 K. If this model were correct, 
then the spectroscopic appearance of the secondary star should show dramatic changes 
over an orbital cycle, oscillating between a state where CO absorption is strong, as in the 
L-dwarfs, to one where methane might dominate, as seen in T-dwarfs. We have obtained 
new phase-resolved H- and K-hand spectra of EF Eri using NIRSPEC on KECK 11, and 
K-hand spectra using NIRI on Gemini-North, to test the irradiated brown dwarf model. 



2. Observations 

EF Eri was observed continuously from 7:21 until 9:00 UT on 24 December 2002 using 
NIRI^ on Gemini- North in queue scheduled mode. For these observations the F/6 camera, 
a four pixel slit and the K-grism were used to produce K-hand spectra with a resolution of 
R ~ 780. Individual spectra, with exposure times of 120 s, were obtained at five different 
positions along the sht. After each exposure, the telescope was offset by 6" along the sht, 
and a new exposure started. In this way forty two individual spectra were obtained in 100 
minutes, covering 123% of an orbital cycle for EF Eri. In addition to this set of continuous 
spectra, four additional spectra were obtained in the period between 6:10 and 6:19 UT. 

To remove telluric features in the spectra of EF Eri, the G2V star HD16358 was 
observed using an identical procedure, but with exposure times of 5 s. At the time of the 
observation, HD16358 had an airmass of 1.44, which is very close to the mean airmass 
attained by EF Eri during the 100 minute observation period. The resulting dataset was 
reduced using a procedure common for reducing infrared spectral data in which the images 
are grouped into sets of five individual exposures. A median of four of these exposures is 

^http://www.us-gemini. noao.edu/sciops/instruments/niri/NIRIIndex.html 
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subtracted from the spectrum of interest to remove the sky contamination. A spectrum 
is then extracted from this sky-subtracted image using the normal longsht spectroscopic 
reduction packages within IRAF, and the process is repeated. Each individual spectrum 
was wavelength calibrated using the spectrum of an argon lamp obtained shortly before our 
observations began. The wavelength calibrated spectra of EF Eri were then divided by the 
median of the five wavelength calibrated spectra of HD16358 to remove telluric absorption 
features. The latter procedure introduces spurious features into the spectra of interest due 
to weak atomic and molecular absorption lines in the G2V spectrum. To remove these weak 
features, we have employed the procedure described by Maiohno, Rieke, & Rieke (1996). 

With only a single standard to work with, the telluric correction for the entire data 
set is not as good as can be produced (see Harrison et al. 2004) when a large number of 
such spectra are available. This is evidenced by dips at 2.0 /im and 2.05 /im in the NIRI 
spectra presented below. The flux calibration of these data was fairly good, and our spectra 
revealed the same large amplitude {/S.K — 0.8 mag) flux variations as reported by Harrison 
et al. (2003). We have assumed that EF Eri has not changed its photometric behavior 
between the time of the last photometric data set (see below) and the epoch of these 
observations, allowing us to use the if-band light curve to produce flux-calibrated spectra. 

EF Eri was observed using NIRSPEC^ on Keck II on 6 September 2003. We used 
NIRSPEC in the low-resolution mode with the 0.38" sht. We employed the four-nod script 
to obtain spectra at four different positions along the sht. For the //-band spectra the 
grating tilt was set to cover the wavelength range 1.54 to 1.83 jim. with a dispersion of 2.91 
A/pix. The if-band setup covered the wavelength range 2.00 to 2.42 //m, with a dispersion 
of 4.27 A/pix. Due to the strong curvature of the spectra produced by NIRSPEC, each 
spectrum has a slightly different wavelength coverage that depends on the position of 



^http: / / www2.keck.hawaii.edu/inst / nirspec / nirspec.html 
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the object along the slit. This effect will be apparent in some of the NIRSPEC spectral 
sequences presented below. The iJ-band data were obtained from 12:29 to 13:57 UT, 
covering 109% of an orbital cycle. The K-hand spectra were obtained in the interval 14:24 to 
15:25 UT, and covered 88% of an orbital period. HD16358 was used as the telluric standard 
for both bandpasses. While the conditions were photometric, we found that the nodding 
routine did not reliably place our targets exactly on the slit during the four-position nod, 
and thus we were unable to extract photometric information from these data. We have 
decided to use the infrared light curves to produce flux-calibrated spectra. 

Three months before the Gemini observations, on 2002 September 24, we were able to 
obtain additional, simultaneous Jii/X photometry of EF Eri using SQIID on the KPNO 2.1 
m. These data were reduced in the same fashion as our earlier observations (see Harrison et 
al. 2003 for full details). The conditions during this run were poorer than those of our 2001 
run, and the J-band light curve (see Fig. 1) has a lower S/N than the earher effort. The 
H- and if-band data, however, reveal identical light curves to the 2001 data set, showing no 
evidence for changes in either the absolute flux levels, or in the amplitude of the variations. 

Harrison et al. (2003) found that the phasing of EF Eri by Bailey et al. (1982), where 
phase corresponded to the appearance of a sharp dip in the J-band light curve, appeared 
to be identical with binary orbital phase. We have used the same set of ephemerides to 
phase the data here. Because of EF Eri's extremely short period, the spectral smearing 
caused by the orbital motion is signiflcant. We have used the evolution code predictions by 
Howell, Nelson, & Rappaport (2001) to estimate the masses of the two components in the 
EF Eri system as Mi = 0.6 Mq and M2 = 0.06 M0. Using an orbital inchnation of 45° (see 
Harrison et al. 2003, and references therein), we estimate the radial velocity semi-amplitude 
of the secondary in EF Eri to be K2 = 340 km s^^. The full orbital variation corresponds to 
only about two resolution elements in the NIRI spectra, but six resolution elements for the 
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Keck spectra. Thus, to allow us to coadd the individual spectra to search for features from 
the secondary star, we have Doppler corrected our data using this estimated value for K2. 

In our ongoing program of continued, occasional monitoring EF Eri, we have obtained 
a new optical spectrum using the EMMI^ spectrograph in the Red Imaging Low Dispersion 
(RILD) mode on the NTT. This spectrum, shown in Figure 2, was obtained on 10 August 
2002, is a 600 s exposure, and has a resolution of 2.3 A. It is interesting to compare this 
spectrum with those presented in Harrison et al. (2003), and Beuermann et al. (2000): 
the Zeeman absorption features are now more pronounced, and the emission from H I has 
continued to decline to the point where it is no longer detectable. 

3. Results 

In the preceding section we noted that we have produced flux-calibrated spectra by 
using the infrared light curves for EF Eri. We did this so as to produce a consistent set of 
data that would enable us to ignore slit losses in either set of spectra. Comparison of data 
sets obtained on widely separated epochs like those analyzed here could be compromised by 
intrinsic variations in EF Eri. While we are confident about the photometric "stabihty" of 
EF Eri up to, and including the Gemini data, there is no assurance that EF Eri was in the 
same state when we observed it with Keck. But to make sense of the Gemini NIRI data has 
required us to assume that the if-band spectra obtained with NIRSPEC represents a similar 
state in EF Eri. The fact that the NIRSPEC if-band data had similar phase-dependent 
morphology gives us confidence in this assumption. 

In Fig. 3 we present the H- and if-band spectra for EF Eri at a number of orbital 
phases. Each of these spectra is the median of four (i/) or five {K) individual spectra. 

^ht t p : // www . Is . eso . org /lasilla/ Telescop es / NE WN TT / emmi / 
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The if-band spectra in this plot have been produced as described above. The telluric 
correction scheme of Maiolino et al. (1996) truncates the blue end of the K-hand at 1.94 
/im. Due to the blueward rise seen in the NIRI spectra of EF Eri, we have also reduced 
the K-hand spectra by simple raw division of these data by the telluric standard, and then 
multiplied the divided spectra by the spectrum of the appropriate blackbody. This allowed 
us to extend our wavelength coverage down to ^ 1.92 fim. This latter method, however, 
introduces weak "emission" features from the G-star division (see Fig.l in Harrison et al. 
2004). Given the low S/N of our data set, the only obvious such feature is the false H I 
Br-7 emission hne at 2.16 /xm visible in the spectra with orbital phases near 0.5. 

3.1. The i/-band Spectra 

During one orbital cycle, the if-band spectra undergo remarkable changes. Near phase 
0.0 (i.e., minimum light), the spectra show a pronounced blueward rise. Note that the data 
have relatively high S/N in this region, as it is not near the blue cutoff of the if-band. 
Alternatively, one could suggest that there is a single, broad absorption feature located at 
1.574 /im. After this dip, the spectrum rises to a broad maximum whose peak is near 1.72 
//m. As we show below, this is not the expected behavior for the spectrum of a normal 
late-type/brown dwarf star. 

As the orbital phase progresses from phase 0.0, the blueward rise rapidly disappears, 
and the flux peak in the spectrum moves to the blue. It is unclear whether this blue feature 
completely disappears, or is simply hidden by the much higher flux levels associated with 
the phase 0.5 cyclotron emission. During the transition to maximum light, a weak rise at 
the extreme red end of the if-band appears. While the S/N of the data in this region is 
quite low, the combination of the H and K-hand data sets suggest that this rise is real, 
and is due to a broad feature that peaks somewhere between the edges of the H and K 
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bandpasses, in the middle of the region strongly affected by telluric water vapor absorption. 
As the orbital phase continues to increase on its way back to phase 0.0, the spectra retrace 
their pre-maximum evolution. 

3.2. The K-hand Spectra 

The if-band spectra obtained with NIRI are shown in Fig. 3, and the K-hand spectra 
obtained with NIRSPEC are presented in Fig. 4. Both data sets show similar evolution 
over an orbital cycle, suggesting that little has changed in the intervening nine months. 
Near phase 0.0, the spectrum shows a continuum that is consistent with that of a late- type 
star, with declining fluxes at both blue and red ends, presumably caused by water vapor 
absorption in the photosphere of the secondary star. The co-addition of the NIRSPEC 
phase 0.0 spectra (Fig. 5) do not confirm this, however, with no clear sign of the expected 
spectral features for a late-type star. This is true even when we correct the spectra for the 
predicted radial velocity of the secondary star. While a feature consistent with ''^^C0(2,o) (at 
2.294 yum) appears to be present, other features with a similar depth do not correspond to 
any of the expected atomic features (e.g., Na I at 2.20 /xm) for a cool star. By coadding 
the entire NIRSPEC data set, the first overtone feature of ^^C0(2,o) appears to remain 
marginally detectable. This would be the first detection of the secondary star in EF Eri. 

As the orbital phase advances, the spectra show rapid changes with the development 
of a pronounced hump blueward of 2.10 /xm, and an increasingly steeper continuum beyond 
2.2 /xm. This blue hump is better defined in the NIRI spectra due to their extended blue 
coverage. Note that the telluric water vapor absorption below 2.0 /xm is very strong, and 
thus the shape of the continuum here is highly sensitive to the telluric correction. While the 
blueward rise begins well before this difficulty, the exact form of the slope of the continuum 
below 2.0 /xm is not well-defined with the current data set. 
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3.3. The Secondary Star of EF Eri 

As we have just noted (see Fig. 5), combining the entire NIRSPEC data set has allowed 
us to tentatively detect the secondary star in EF Eri. While the S/N of these data are not 
high, if the secondary star were relatively normal, it should be more prominent than seen 
here. By modeling the optical spectrum of EF Eri, Beuermann et al. (2000) have derived 
distances for EF Eri in the range 80 < d < 128 pc, while Thorstensen (2003) has measured 
a distance of 163 ± 58 pc via parallax. Thorstensen finds that an astrometric solution 
with a distance of 113 pc is possible by adjusting some of the input criteria. As shown in 
Harrison et al. (2003), the optical photometry of EF Eri is well-modeled by a 9500 K white 
dwarf, and this white dwarf provides no more than 10% of the flux in the if-band. If we 
assume that the secondary star is responsible for the remainder of the K-hernd flux at phase 
0.0 {K = 15.55), the absolute magnitude of EF Eri is in the range 10.0 < Mk < 11.0. This 
absolute magnitude corresponds to a spectral type in the range M8 to LI (ignoring any 
effects due to the Roche lobe filling nature of this object). 

A secondary star with a spectral type near M9 should have strong CO features, as well 
as a detectable Na 1 doublet. We have compared the combined H- and K-hand spectra of 
EF Eri with published spectra for stars and brown dwarfs in the range M7V to T8. While 
the overall slope of the continuum of the phase 0.0 K-hand spectrum is consistent with a 
very late- type M dwarf (see Fig. 6), the //-band spectrum bears no resemblence to such an 
object. The departure becomes more pronounced when we compare the spectrum of EF 
Eri to those of L-dwarfs. While brown dwarfs around L5 show a similar spectral slope at 
1.55— 1.6 //m, there is no feature comparable to the strong flux peak at 1.7 /um. While the 
absence of Na 1 in the K-hand is consistent with a spectral type later than LO (McLean et 
al. 2003), the water vapor absorption bands of such an object are much deeper and depress 
the blue side of the K-hand spectrum much further to the red (> 2.1 /im) than is seen in 
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EF Eri. Also, the first overtone CO bandhead at 2.29 //m remains strong in brown dwarfs 
until CH4 takes over in the late L types. 

In addition to the difference in H- and X-band spectral features, the JHK photometry 
is almost impossible to reconcile with any known dwarf. While the H to i^^band flux ratios 
for EF Eri are more-or-less consistent with those of late type objects, giving rise to [H — K) 
— 0.42 (ps M8V), the observed (J — K) color (> 2.5) exceeds even the reddest IR colors 
observed in mid-L dwarfs (Legget et al. 2002). 

3.3.1. Model Spectra for Carhon/ Oxygen- depleted Atmospheres 

As shown in Harrison et al. (2004), many CV secondary stars show evidence for 
weak CO absorption features. Since the water vapor features in those objects appear to 
be relatively normal, Harrison et al. suspect that Carbon is deficient in those systems. 
This view is supported by analysis of FUV observations of CVs by Gansicke et al. (2003) 
where large N V/C IV ratios have been found, suggesting Nitrogen enhancements, and 
Carbon deficits in the photospheres of the white dwarf primaries, presumably arising from 
matter transferred from the secondary star. It would therefore not be unexpected for the 
secondary star in EF Eri to have anomalously weak CO absorption features. A lowered 
Carbon abundance would also reduce the CH4 absorption expected for late-L and T dwarfs, 
and thus could produce much redder (J — K) colors through the range of possible spectral 
types for the secondary of EF Eri. 

To test this hypothesis we have calculated theoretical spectra using the general stellar 
atmosphere code PHOENIX (Hauschildt et al. 1999) for a variety of chemical compositions 

that might be expected in a CNO-proccsscd stellar core. Depth-dependent partial pressures 
for more than 650 species are calculated using the EOS of AUard & Hauschildt (1995) 
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by solving chemical equilibrium for 40 elements, as described in AUard et al. (2001), 
including dust formation. Model structures and flux spectra are computed self-consistent ly, 
treating line blanketing by the direct opacity sampling method, with line strengths taken 
from a master database of 43 million atomic and 700 million molecular lines. In addition, 
dust extinction profiles for 31 grain species are included. Major updates relative to the 
AMES-Cond and AMES-Dusty models of Allard et al. (2001) arc the addition of 30 million 
methane lines from Homeier et al. (2004a), and the explicit treatment of the gravitational 
settling of dust grains developed by Allard et al. (2003). In addition, the Ames list of water 
line opacities from Partridge & Schwenke (1997) has been replaced by the earlier list from 
Miller et al. (1994) for these models. The latter list is less complete for high temperature 
steam, but has been found to give better agreement with the IR spectra of cool brown 
dwarfs (Homeier et al. 2004b). Nevertheless we emphasize that none of the currently 
available databases for molecular opacities accurately reproduces the observed shape of all 
absorption bands in the IR. These models also suffer from incomplete or missing data for 
several molecules, such as CrH and FeH. For example, the ff-band spectra of mid-L dwarfs 
show a moderately strong absorption feature near 1.6 ^m (see Fig. 6) that is not reproduced 
by any of our models. Several features in this region have been identified by Gushing et al. 
(2003) with the E^U-A^U system of FeH, for which no hne opacity data are available at 
this time. Possibly due to this missing opacity source, and other uncertainties in input data 
such as H2O opacities and dust extinction profiles, our theoretical spectra tend to produce 
{H — K) colors that are too blue when compared to observed L dwarfs (Legget et al. 2001). 

We have constructed grids of models based on revised solar abundances as compiled by 
Lodders (2003), except for Carbon, Oxygen and Nitrogen, which were altered as indicated 
for the individual models. The basic atmospheric parameters span the range 1400 K < 
Teff < 2500 K, and \ogg= (4.5, 5.0, 5.5). First, investigating the effects of lowering the 
Carbon abundance to [C/Fe] = (— 0.5, —1.0, —1.5), we find that a depletion to at least 
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0.1 X solar is required to explain the weakness or absence of the CO overtone bands. Prom 
synthetic photometry we can also exclude models with T^q significantly higher than 2000 K, 
as these disagree with the observed J magnitudes even assuming the upper limit of 5.5 for 
the distance modulus. Below 1600 K, the strong methane bands at 2.2- 2.3 //m become 
prominent even at [C/Fe] = —1.0, but lowering the Carbon abundance further does not 
improve the fit to the observed i^-band spectrum at such low temperatures. We find best 
agreement for the red part of the K-hand spectrum for models with 1600 < Teff < 1800 K, 
as shown by the bottom spectrum in Fig. 7. But a Carbon deficiency alone cannot explain 
the early onset of absorption to the blue side of the K-band, which appears at much shorter 
wavelengths than typically seen in both our models and in observations of L dwarfs. Nor 
can a lack of Carbon be responsible for the peculiar shape of the ff-band spectra. The 
main opacity source at these wavelengths is water vapor, and significant changes can only 
be expected if the Oxygen abundance is reduced as well. Such effects might be expected if 
the secondary CNO cycle has been active for some time, though models by Marks & Sarna 
(1998) do not indicate significant changes in the Oxygen abundance for CNO processed 
material. We have nonetheless produced test models for [0/Fe] = —0.5 to explore this 
case. As one would expect, these models (Fig. 7) show an onset of the water vapor bands 
both in the H and K-heaid at shorter wavelengths, with a K-hand spectrum that fits our 
observations quite well. For comparison, a model with additionally enhanced Nitrogen 
abundance, as would be expected with the depleted Carbon and Oxygen transformed into 
^^N, is also shown. 

Synthetic colors calculated for these models show that the Nitrogen-enriched 
atmosphere produces the closest match to the {J—K) and {H—K) photometry. We have not 
attempted to modify our input chemistry further to optimize these fits, since such results 
would necessarily be of limited significance given the uncertainties in the available opacity 
data and the limited wavelength coverage of our spectra. Clearly a direct detection of the 
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water absorption in the iZ-band, which we expect to occur at 1.5- 1.55 yum, and J-band 
spectra, would help in resolving these uncertainties. Such data would also allow for the 
study of irradiation, which is expected to have the strongest impact at shorter wavelengths 
(Barman et al. 2000). In conclusion, a brown dwarf with a considerable deficiency of 
Carbon, a lesser depletion of Oxygen, and an enhancement of Nitrogen, reproduces much 
of the secondary's IR SED, and the shape of the i^-band at phase 0.0. The additional 
absorption at 1.6/Ltm cannot be accounted for in the current models. If this band is indeed 
due to FeH, we expect it to become relatively stronger when the steam bands are reduced. 
In that case the spectra could probably reproduce the "base" of the ff-band emission. The 
shape of the actual peak, however, does not correspond to any known stellar features, and 
also shows changes with phase that cannot be explained by irradiation effects. We conclude 
that cyclotron emission at 1.6-1.8/um contaminates the infrared spectrum for all orbital 
phases. 

3.4. The Cyclotron Features 

It is clear that the H-hsmd spectra are inconsistent with that of a normal (or chemically 
altered) late-type or brown dwarf secondary star. If the phase 0.0 dip near 1.6 /im is 
associated with FeH, then we estimate that 30% of the if-band flux at this juncture is 
due to cyclotron emission. Thus, to produce the observed, phase-dependent changes in 
the H-hand spectra probably requires two separate cyclotron features in this bandpass. If 
we assume the phasing and geometry from Bailey et al. (1982), the phase 0.0 cyclotron 
features that we see in the if-band are presumably associated with the pole that accretes 
during the high state. This suggests that the single, dominant cylcotron feature at phase 
0.5 might be from the other pole. 

While the minimum light K-hand spectrum is fairly consistent with the shape of a cool 
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star spectrum, with little evidence for cyclotron emission, the maximum light spectrum is 
clearly contaminated by at least one, possibly two cyclotron features. Obviously, there is a 
cyclotron harmonic that peaks below 2.0 /im. The rapid change in slope of the red half of 
the K-hand spectrum, and the dip near 2.07 yum, suggest the possibility of an additional, 
weaker cyclotron feature at phase 0.5. To better examine the phase 0.5 cyclotron features, 
we have subtracted the phase 0.0 spectrum from the phase 0.50 spectrum, to produce 
the residual spectrum shown in Fig. 8. The result shows two strong cyclotron features, 
one that peaks in the if-band near 1.58 /im, and one that (we estimate) peaks near 1.93 
/xm. The possibility of a second cyclotron feature in the if-band seems more likely after 
this subtraction. If the secondary star has a similar spectrum at phases 0.0 and 0.5, then 
any water vapor feature should disappear upon subtraction. Ignoring the false H I Br-7 
emission feature, there continues to be a modest dip near 2.07 fj,m and a small rise that 
peaks near 2.15 /xm and that falls off slowly beyond 2.2 //m. 



4. Discussion 



Regardless of the exact interpretation of the spectra presented here, it is clear that the 
model proposed by Harrison et al. (2003) is incorrect. The photometric modulations in the 
if-band are certainly being driven by the changing strength of the cyclotron feature at the 
blue end of this bandpass (and possibly aided by a broader, and weaker feature that peaks 
near 2.15 //m). There is very little evidence for the secondary star in our K-hand spectra, 
with only the suggestion of a weak first overtone CO feature in the median of the if-band 
data set. While the presence of weak CO features in a CV is not surprising (see Harrison 
et al. 2004), stars with the late spectral type expected for the secondary in EF Eri have 
if-band spectra that are completely different to what we have observed. While model brown 
dwarf spectra with dramatic alterations in the abundances of the major opacity sources at 
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these wavelengths do a somewhat better job at fitting the data, neither the observations or 
the models are currently adequate enough to truly investigate the nature of the underlying 
secondary star. 

The fact that cyclotron emission can still dominate the spectrum of EF Eri in this low 
state is surprising given that EF Eri exhibits no outward signs of continued accretion. The 
optical and infrared spectra lack Balmer and Brackett emission from H I. The optical BVRI 
spectral energy distribution and spectrum are completely consistent with that of an inactive 
magnetic white dwarf. The model white dwarf atmosphere, with Tcfr = 9500 K, used by 
Harrison et al. (2003) to compare with their 2001 March 11 spectrum of EF Eri, does an 
excellent job at fitting the 2002 August 10 optical spectrum we obtained with the NTT. 
Clearly, the white dwarf photosphere alone is too cool to supply significant quantities of 
ionizing photons. It appears that the white dwarfs in ultra-short period, low mass transfer 
rate CVs are significantly cooler than those in longer period, or higher mass transfer rate 
systems (c.f., Szkody et al. 2002). 

Four magnetic CV systems have now been found that display a single, dominant 
cyclotron feature in their optical spectra that have been attributed to very low mass 
transfer rates: SDSS J155331. 12+551614.5 and SDSS J132411. 57+032050. 5 (Szkody ct 
al. 2003), RBS0206 (= CQ Get; Schwope, Schwarz, & Greiner 1999), and HS 1023+3900 
(= WX LMi; Reimers, Hagen & Hopp 1999). Along with the single very large cyclotron 
feature is evidence for a second, weaker feature at the next lower /higher harmonic number 
covered by their optical spectra. Of these four polars, both SDSS J155331. 12+551614. 5 
and J132411. 57+032050. 5 have optical hght curves that resemble the infrared hght curves 
of EF Eri (see Figs. 2 and 6 of Szkody et al. 2003). SDSS1553 shows optical variations 
with peak-to-peak amplitudes of Am ^ 1 mag in Fand i?, while SDSS1324 has variations 
in the Sloan r-band of Ar =1.3 mags. Both sets of light curves are roughly sinusoidal. 
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like the H and if-band light curves of EF Eri. Unhke EF Eri, however, the spectra of both 
SDSS1553 (Fig. 1 in Szkody et al. 2003) and SDSS1324 (G. Schmidt, priv. comm.) reveal 
H I emission. 

The interpretation of the spectra of SDSS1553, SDSS1324, CQ Get and WX LMi is 
that they have entered phases of extremely low accretion rates, the so-called "bombardment 
solution" (Kuijpers & Pringle 1982), where any heating is from particle collisions instead of 
shocks. This state is characterized by low plasma temperatures, and mass transfer rates of 
M < 10^^^ M0 yr~^. In each of these four polars, the harmonic numbers of the dominant 
cyclotron features are in the range 2 < n < 4, with magnetic field strengths above 45 MG. 
Magnetic field strength estimates for EF Eri range from 10 to 21 MG. Assuming B ~ 14 
MG, the n = 4 harmonic would peak at 1.96 /xm, and the n = 5 harmonic would peak at 
1.58 /xm. The n =3 harmonic would peak near 2.6 /xm, while the n = 6 and 7 harmonics 
would appear in the J-band. It is worth noting that during the H- and if-band maxima, 
the J-band light curve shows a minimum. Given this fact, and the need for very high 
accretion rates to produce power in these higher harmonics, it is likely that the n = 6 and 7 
harmonics are not significant sources of flux in EF Eri. The fact that the J-band light curve 
has a maximum during the H- and K-hand minima suggests activity from another pole with 
a signiflcantly higher magnetic field strength. The flux peak at 1.72 ^m, and the upturn at 
the blue end of the phase 0.0 iJ-band spectra, along with the lack of significant cyclotron 
emission in the JT-band indicate that the phase 0.0 cyclotron emission comes from a region 
with a significantly different field strength. A fit to the cyclotron features of the high state 
JHK spectra of EF Eri by Ferrario ct al. (1996) suggested emission from two regions with 
magnetic field strengths of 16.5 and 21 MG. An accreting pole with a field strength of 21 
MG, would not produce a cyclotron feature in the X-band, but would produce harmonics 
at 1.7 and 1.3 //m, consistent with our phase 0.0 spectra. 
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That the magnetic field structure of EF Eri is complex, including strong evidence for 
emission from more than one pole, has been demonstrated by Hutchings et al. (1982), 
Piirola et al. (1987), Meggit & Wickramasinghe (1989), and Ferrario et al. (1996). Recent 
"magnetic field topology" mapping of the magnetic field of the white dwarf in EF Eri 
by Reinsch et al. (2003) reveals this complexity: there are two main polar caps with 
field strengths as high 100 MG at their centers, as well as evidence for multiple, zonal 
components. The complex structure of the field in EF Eri must be partly to blame for the 
results presented here. Modeling the cyclotron features that would arise from accretion 
onto such a field would be helpful in unraveling our spectroscopic results. 

The lack of detectable H I emission from EF Eri, when compared to the other members 
of its small family, is interesting. It is highly likely that the current accretion rate in EF 
Eri is lower than the other four objects, and is simply not sufficient to supply the heating 
necessary to generate the UV photons required to create H I emission. In the bombardment 
scenario, the accretion shock disappears, and the white dwarf photosphere is directly heated 
by particle collisions. A lower accretion rate means less heating, and therefore less ionizing 
radiation. 

The secondary star of EF Eri remains elusive. Beucrmann et al. (2000) have presented 
arguments that the spectral type of the secondary star in EF Eri must be later than 
M9, with the most probable spectral type later than L3. We can estimate the maximum 
amount of fiux from the secondary by assuming that the fiux level of the dip seen in the 
phase 0.0 H-hand spectrum (at 1.57 /xm) represents the continuum level of the underlying 
secondary star. This corresponds to /4ec > 16.5. Assuming 80 < d < 128 pc, the absolute 
if magnitude of the secondary star is in the range 11.0 < Mh < 12.0, consistent with an 
early/mid-L dwarf. Thus, EF Eri joins WZ Sge (see Howell et al. 2004) as the cataclysmic 
variables with the strongest evidence for harboring brown dwarf-like secondary stars. 



5. Conclusions 



We have obtained phase-resolved infrared spectroscopy of EF Eri and conclude that 
the large amplitude variations seen in its infrared light curve are not due to a heated brown 
dwarf, but has its origin in cyclotron emission. EF Eri appears to join the small family of 
polars where the accretion rate has dropped to a very low level. Unfortunately, our data are 
inadequate to allow us unravel the complex cyclotron emission that is present in EF Eri. 
Phase-resolved J-band spectroscopy is required to make further progress on disentangling 
the cyclotron spectrum of EF Eri. Given that the mean J-band magnitude of EF Eri is 
17.3, multiple orbits of data on 8 m-class telescopes will be needed to obtain spectra with 
sufficient S/N. To fully interpret the spectroscopic observations, phase- resolved, broad-band 
JHK polarimetric observations of EF Eri will be essential. It would also be extremely useful 
to have a radial velocity curve for the EF Eri system to allow us to properly phase, and 
coadd the infrared spectra to search for features from the secondary star. Such data could 
be acquired through optical spectroscopy of the Zeeman absorption features of H I. 
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Fig. 1. — The infrared light curves of EF Eri obtained with SQIID on the KPNO 2.1 m. 
The absolute fluxes and light curve amplitudes have not changed when compared to those 
presented in Harrison et al. (2003). The data sets have been plotted over two orbital cycles 
for clarity. 
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Fig. 2 — An optical spectrum of EF Eri obtained with EMMI on the NTT on 2002 August 



10. The accretion rate in EF Eri has dechned to the point where there is no longer detectable 
emission at Hct. The features from the Zeeman splitting of the H I lines, however, are 



clearly seen. The central, unshifted absorption components of H I are weaker than expected. 



suggesting the possibility of continued low-level emission. 
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Fig. 3.— The Keck NIRSPEC i/-band, and Gemini NIRI K-hand, spectra of EF Eri. Both 
data sets have been flux-cahbrated using the infrared hght curves. For clarity, the H- and 
if-band spectra for each orbital phase have been offset in flux. The mean orbital phases of 
the spectra are hsted. 
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Fig. 4. — The NIRSPEC if-band phase-resolved spectra of EF Eri. These data have been 
flux-calibrated using the infrared hght curves and then offset for clarity. The same phase- 
dependent morphology seen in the Gemini data set is present here. Each spectrum is the 
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Fig. 5. — The median of the entire NIRSPEC if-band data set (top, offset in flux), and a 
smoothed (to FWHM = 60 A) median of the seven spectra closest to phase 0.0 (bottom). 
A possible, very weak absorption feature at the position of the first overtone of ^^C0(2,o) (at 
2.294 /im) is seen in both spectra, and appears to be the only indication of the presence of 
the secondary star, besides the decline at the red end of the if-band that might be associated 
with water vapor absorption. 
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Fig. 6. — Comparison of the phase 0.0 spectrum of EF Eri to those of M9 and L5 dwarfs 
(from McLean et al. 2003). The sharp feature at 2.294 /i,m in the comparison spectra is due 
to the first overtone feature of CO. 
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Fig. 7. — Comparison of the phase 0.0 spectrum of EF Eri to brown dwarf model spectra 
generated using PHOENIX. The temperature of each model is listed, and each of these models 
has a gravity of logg = 5.0. Below the temperature are the abundances of Carbon, Nitrogen, 
and Oxygen (in dex, vs. solar). The top model has deficits of Carbon and Oxygen, and 
enhanced levels of Nitrogen (4x solar). The bottom- most model simply has a Carbon deficit 
(O.lx solar). The Oxygen deficient models fit the observed iC-band spectrum due to their 
shallower water vapor absorption. None of the models we developed, however, could explain 
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Fig. 8. — The residual spectrum created by subtracting the phase 0.0 spectra from the phase 
0.5 spectra. The result is two apparent cyclotron humps that peak near 1.58 and 1.93 /im. 



